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Abstract. To improve convergence results obtained using a framework for 
unsymmetric meshless methods due to Scliaback (Preprint Gottingen 2006), 
we extend, in two directions, the Sobolev bound due to Arcangeli et ah (Numer 
Math 107, 181-211, 2007), which itself extends two others due to Wendland and 
Rieger (Numer Math 101, 643-662, 2005) and Madych (J. Approx Theory 142, 
116-128, 2006). The first is to incorporate discrete samples of arbitrary order 
derivatives into the bound, which are used to obtain higher order convergence 
in higher order Sobolev norms. The second is to optimally bound fractional 
order Sobolev semi-norms, which are used to obtain more optimal convergence 
rates when solving problems requiring fractional order Sobolev spaces, notably 
inhomogeneous boundary value problems. 



1. Introduction 

Over the past few years, increasingly general bounds of Sobolev semi-norms, in 
terms of discrete samples, have appeared. Such bounds are often called sampling 
inequalities. A rather general sampling inequality was established by Arcangeli et 
al. [T] , and is stated as Theorem II. 1[ with notation given in Section 11.11 

Theorem 1.1. [1, Theorem 4.1] Let H. be a Lipschitz domain in R", so that the 
domain D, satisfies the cone property [U Page 185] with radius p > and angle 
9 e (0,7r/2]. Furthermore, let p,q,>c G [l,c)o] and let r he a real number such 
that r > n, if p ^ I, r > n/p if 1 < p < oo, or r G N*, if p = oo. Let 
Iq = r — n{l/p— ^/q)^ and 7 = max{p,q,>f}. Then, there exist two positive 
constants X)r {dependent on 9,p,n and r) and C {dependent on fl,n,r,p,q and >!r) 
satisfying the following property: for any set A d fl (or A C if p = 1 and r — n) 
such that d — 6 (A, ft) < dr ( c.f. ([2]) ) for any u £ W^''^ {Vi) and for any real number 
I satisfying 1 = 0,..., Imax, we have 

(1) < C (d-'-(Vp-i/.)+ 1^1^^^^^ + d-h-^ II u\,\\^) . 

where /max '-^ \h)\ — 1, unless the following additional conditions hold, in which 
case Imax Iq-' r G N* and either (i) p < q < 00 and Iq G N, (ii) {p,q) = (l,oo), 
or (Hi) p > q. 

This sampling inequality generalizes those of Madych [6 and Wendland and 
Rieger |15) . by greatly extending the range of parameters r,p,l, and >f. While 
Theorem 11.11 applies to functions with finite smoothness, an analogous bound for 
functions with infinite smoothness has been provided by Rieger and Zwicknagl [10] 
which achieves exponential factors. 

Arcangeli, et al. [l] used Theorem 11.11 to derive error bounds for interpolating 
and smoothing (m, s)-splines, an application which we do not consider. Instead we 

1 



2 



ANDREW CORRIGAN, JOHN WALLIN, AND THOMAS WANNER 



are interested in another major application of these Sobolev estimates: Schaback's 
framework for unsymmetric meshless methods for operator equations [12] , see also 
the earlier version '13] . A sampling inequality is necessary for unsymmetric meshless 
methods, such as Schaback's modification of Kansa's method J^ilSj, which involve 
an overdetermined system of equations in general. In an attempt to improve the 
order of convergence obtained using Schaback's framework, we extend the bound 
of Arcangeli, et al. in two ways. 

Our first extension is to loosen the restriction ^ G N to allow for fractional order 
Sobolev norms on the left hand side of the sampling inequality. In the context 
of Schaback's framework, this will result in more optimal convergence results in 
terms of both the test and trial discretization parameters. Otherwise, the test 
discretization would require a higher rate of refinement. 

Our second extension is to incorporate discrete samples of arbitrary order deriva- 
tives into the bound. The reason for this is that ([T|) has a factor d^' in its second 
term which is insufficient for achieving a uniformly stable test discretization for 
higher order Sobolev norms in Schaback's framework. With this modification to 
incorporate samples of higher order derivatives we will be able to come closer to 
achieving such a test discretization, resulting in higher order convergence results. 
This introduces a new parameter /i, for which previous sampling inequalities coin- 
cide with the choice /i = 0. 

1.1. Notation. We employ the notation of Arcangeli, et al. [U Section 2]. 

As in [H Section 4], we assume throughout this paper that f2 is a bounded 
domain in M" with a Lipschitz-continuous boundary, so that the domain satisfies 
the cone property JL, Page 185] with radius p > and angle 6 E (0,7r/2]. For a 
given finite subset A of f2, the fill distance is defined as 



The following restates a portion of their notation. For all r G [0, oo] and p G 
[l,oo], the Sobolev norm is denoted by |H|,. p^, while the Sobolev semi- norm is 
denoted by H^p^^. The set N* = {1, 2, 3, ...},' while N = {0,1,2,...}. The space of 
polynomials over R" with degree less than or equal to k is denoted by Pk- 

We make the following additions to their notation. Let (il) :— W^''^ (il) and 



Given a function v G W^''' (il), the vector- valued function consisting of its par- 
tial derivatives is denoted by Dv. The surface area of the n-dimensional ball is 
denoted by IS"""^]. In Section O a generic constant C appears in many proofs, 
whose particular value may change, but with the parameters on which it depends 
either indicated in parentheses or stated explicitly in the exposition. We will often 
substitute dependencies with others, possibly taking the maximum or minimum 
value, as required by its application, of the constant over a finite range of values. 
In Section [31 only the dependence of constants on the discretization parameters r 
and s is explicitly stated since we regard the spaces and mappings in that section 
as fixed. 



(2) 





A SAMPLING INEQUALITY 



3 



2. Extension of the Sobolev Bound 



2.1. Fractional Order Sobolev Spaces. This section concerns fractional order 
Sobolev norms and the results of this section will be used to generalize [U Proposi- 
tion 3.4] to Proposition l2.7l Lemmas 12.21 and [ 2 .31 each require an extension operator 
which satisfies ^ for zero-average functions over a ball. This property is not pro- 
vided by standard extension operators since they involve a domain-dependent con- 
stant and the full Sobolev norm in the bound, rather than a domain-independent 
constant and the Sobolev semi-norm. 

Lemma 2.1. If q (z [1, oo], r > 0, and xq S M" then there exists a linear, continu- 
ous operator 

E : W^^'^ {B{xo,r)) ^ W^^" (M") 
such that for all v G W^^'* {B {xo,r)) 
(3) Ev = V a.e. in B (xq, r) , 

and 
(4) 



// in addition v e (^B{xo,r)^ then Ev £ (R"). 

Proof. Let v G W^'"^ {B (xq, r)) and v := v o F where F : x ^ rx + x^j. From a 
change of variables it follows that v G W"^'^ {B (0, 1)) with semi-norm 

From 5j Section 5.4, Theorem 1], there exists a linear, continuous extension oper- 
ator 

E-.W^^'^ (B(0,l))^M^i'« (K") 
such that for each v G W^''^ {B (0, 1)) 



and 



Ev 



Ev = V a.e. in B (0, 1) , 



where the dependence on n is through B (0, 1). That 

l!^lk9,s(o,i) - C'('^'9)l*k9.s(oa)' 

follows from specializing a Poincare inequality given in [5, Section 5.8, Theorem 2] 
to the unit ball and that 



(5) 

Let 
(6) 



s(oa) 



w = 0. 



B{xo,r) 



Ev := 



(ev) 



so that the result ^ holds. From another change of variables, it follows that 



\Ev\ 



1,9,1 



Ev 



1,9,1 
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The result ([4]) follows by combining the preceding relations. Finally, from the 
proof of O Section 5.4, Theorem 1] it follows that if {) e (^B (0, 1)^ , then 

Ev e (M"), so that if V e (s(a;o,r)) then Ev e (R"). □ 

Based on this extension, we obtain Lemma 12. 2[ which is similar to a result used 
by Bourgain et al. '2, Eq. 2], but uses a domain- independent constant and semi- 
norm in the bound. 

Lemma 2.2. If q e [l,oo), h e M", and v e W'^-'i {B {xo,r)) then 
\Ev {x + h) - Ev (x)]" dx] <C (n, q) \h\ 

lB{xo,r) J 

Proof. First suppose that v is in the subset 



(7) 



C 



1 (^B{xo,r))nW''^B{xo,r)) 



which is dense in W^''^ {B (a;o,r)). From Lemma [2?11 it follows that 



\Ev [x + h) - Ev {x)\'^ dx < 



B{xo,r) 



\Ev (x + h) - Ev {x)\''dx 
"1 d 



I I ^Ev{x + th)dt 
JK" Jo 



dx 



< 



dt 



Ev {x + th) 



dtdx 



= [ f \DEv{x + th)-h\'^ dtdx 

JWL" JO 

< \h\'' [ [ \DEv{x + th)\Utdx 
Jr" Jo 

' [ [ \DEv{x + th)\'^ dxdt 
Jo Jr" 



< \h\'' j \Ev{- + th)\l^ ,,^^dt 

= \K\Ev\%,Rr. 

< C{n,q)\h\''\v\l^^^^^^^^^ 

From this, the result follows for all fmictions in W^'"^ {B {xo,r)) via a standard 
density argument. □ 

Lemma 2.3. If x,x + h £ B {xq, r), and v E 14^^'°° [B {xq, r)) then 
\v{x + h)-v{x)\<C{n)\h\\v\,^^,,^^^,.^ 

Proof. In the proof of [SI Section 5.8.2, Theorem 4] it is shown that w is a Lipschitz 
function with constant \Ev\^ ^ jj„, where the extension operator constructed in [3 
Section 5.4, Theorem 1] is used. However, the extension operator from Lemma |2. II 
could be substituted so that the result then follows by ([4]). □ 

In the bound provided by Lemma [2^ the factor r^~'^ will be the key to general- 
izing sampling inequalities to optimally bound fractional order semi-norms. 
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Proposition 2.4. If q e [I, oc], e e (0,1), and v e W^'*^ {B {xo,r)) then 
(8) \<,M.<ur) < C (n, q) (1 - e)-^/^ r^- . 

Proof. Since the semi-norms tliat appear in ([8]) are invariant with respect to a shift 
in value of w by a constant, it suffices to only consider v G W^''' {B [xq, r)). 
Case q G [1, oo): Let y £ B {xo,r) and 

B {xq, r)-y:= {x~y : x G B {xo,r)} , 

so that B (xo, r) - y C S (0, 2r) . 



K,q,BM = / / '7^' '^^^^ 

JB(x„.r)JB{xo,r) \x - v\ 



dh 



B{xa,r) J B(xa,r)~y \n\ 

\Eviy + h)-Ev iy)\\ 

. TT^q dh 

B{xo,r) J B{0,2r) |"| 

SB(xo,r)\E^ {y + h)-Ev {y)\''dy 

'^("•'g) I en-l I (-9r^(l-')« l„l« 
(1 P M ^ l^^ll,g,B(xo,r) • 



< 



< 



Case q = oo: 



. |w (a;) — w (w)| 

x,y^B(xa,r),x^y i/l 



< 2C{n,q)r'-^\v\,^^^g^^^^^y 



□ 



Remark 2.5. The explicit constant (1 — e) which blows up as e increases 

towards one, is a manifestation of the "defect" of intrinsic fractional order Sobolev 
semi-norms studied by Bourgain et al. 2\. 

Corollary 2.6. // g e [1, oo], ^ e [0, cx)], and v eW''^'^^'> {B {xo,r)) then 

(9) \^kq,BM < C {n, q, ll\)K{\l]~l, q) r ^'1 \v\^,^ , 
where 

f 1 for / G N or (7 = cx) 

(10) if(R1-^g):=| for;^i^„,rf^<^ ■ 
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2.2. An Auxiliary Result. The following result applies Corollarv l2.6l to general- 
ize [U Proposition 3.4]. 

Proposition 2.7. Let p,q, >c [1, oo] such that p < q. Let r he a real number such 
that r > n/p, if p > 1, or r > n, if p = 1. Finally, let k = [rl — 1, A = dimP^, 
and Iq = r — n/p+ n/q. Then, there exists a constant R> 1 (dependent on n and 
r) and, for any M' > 1, there exists two constants C (dependent on M' ,n,r,p,q, 
and >f) and K > I (explicitly dependent on \l~\ — I and q, of. (I10|) l. satisfying the 
following property: for any d > and any t G R", the open ball B (t, Rd) contains 
K closed balls Bi, . . . Ba of radius d such that, for any v G W^'^ (S (t, M'Rd)) , for 
any b G ^f^iBi and I G [0, Ima^], 

(11) H.,.mM'm) <CK IH,^^ + d"/'-' II «|J| J , 

where we have let /max '■— \lo \ — 1; or Zniax := if the following additional conditions 
hold: r G N* and either (i) p < q < oo and Iq G N, (ii) {p,q) = {l,oo), or (Hi) 
1 < p = q < oo. 

Proof. The case that / G N is established by [ll Proposition 3.4]. Suppose that 
I ^ N. The hypotheses imply that /max G N, so that [/] < /max, and thus the 
result follows by combining (fTTj) for / = [/] with Corollary 12. 6( using the fact that 
(Af'i?)'^'^"' < M'R, and substituting the dependence on [/J and R with n,r,p and 
q. □ 

2.3. Sobolev Bounds. For p < q, the following result generalizes [2 Theorem 
4.1] to bound fractional order Sobolev semi- norms. No generalization to bound 
fractional order Sobolev semi-norms is made for p > q, since we have not obtained 
the relation [U Eq. 2.1] for the case that / is fractional. 

Theorem 2.8. Let p,q, x £ [1, oo] and let r be a real number and fj, a nonnegative 
integer such that r — fi > n, if p = 1, r ~ fi > n/p if 1 < p < oo, or r — fi £ N* 
if p — oo. Let Iq = r ~ fi ~ n{l/p— ^/q)j^ and 7 = max {p, g, . Then, there 
exist three positive constants 'Qr (dependent on 9,p,n,r and /i), C (dependent on 
Q,,n,r,p,q, and k), and K > 1 (explicitly dependent on [/] — / and q, cf. (|10p ). 
satisfying the following property: for any set A C (or A C if p = 1 and 
r — = n) such that d = 5 (A, fj) < 0,,, for any u G W"^'^ {Vl) and, if p < q then for 
any real number I G [0, /max], otherwise if p > q then for any integer I — 0, . . . ,l„ 



^max; 



(12) |u|,_^_^ <c k[ d^-'-"(i/P-i/.)+ 1^1^^^^^ + d"/7+P-; 



n ^""1 



where we have letl^^ia^ '■— f/ol^li or /max '■— /o if the following additional conditions 
hold: r G N* and either (i) p < q < 00 and Iq G N, (ii) {p,q) — (1,cxd), or (Hi) 

p>q- 

Proof. The case that fi — and / G N is Theorem 4.1]. The proof of this theorem 
for fi — and / N can be obtained by reusing the proof of [1] Theorem 4.1], but 
applying Proposition 12 . 71 instead of [TJ Proposition 3.4], which allows for / to be of 
fractional order for p < q, and introduces the constant K. We now consider the case 
fj. > 0. Let a be a multi- index such that \a\ — and therefore d"u G W'^~'^'P (fJ). 
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It follows from the case that /x = that in the situation required by the present 
hypotheses 

We have for all a satisfying |q!| = /i that 



and 



\d"nUL< 



1/31=^ 



The result follows immediately for q — oo. Otherwise, if 1 < 9 < 00, using that 



H — fi.q.^i 



the results follows with an additional factor {a : \a\ — /i})^^^ in the constant C, 
whose dependence on /i can be substituted with n, r and p. □ 



Corollary 2.9. Given the situation of Theorem \2.8\ with a constant dr now depen- 
dent on 9,p,n,r,p and q, and the additional assumption that r — / G N, then we 
have 



n 9" 



l"l<A' 



Proof. From Theorem 12.81 there exists three positive constants i:)r {0, p,n,r,p, q), 
C {Q,n,r,p,q, h), and K — l,q) > 1, cf. (fTO]) . such that for d < dr and r] = 

0, . . . , L^J 



+d 



'n/j+irj+ti-ll]) 



n 



|a| = (»)+M-W) + 



and for 77 = Z that 



n ^"^1 

\a\=fj. 



We have taken the constants to be the minimum or maximum over rj — 0, . . . , [l\ ,1 
as required. Additionally we have restricted Or to be at most one. For rj = 0, . . . , \ 
we have applied Theorem 12 . 81 with r = r — l-\-r] and fj, — {rj + fi — introducing 
a dependence of dr on p and q through /, which along with n and r has substituted 
for the dependence on /i. It also follows for all 77 = 0, ... , [l\ that 



(13) 



n 



go 



|Q| = (r,+M-L'J) + 



< 


n ^""Ia 




lal<M 
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It follows from 



with a similar bound holding for 
and Or < 1 that for all d <()r, 

^j^4) i^"/7+('j+A'-L'J)+-f) < (f-h+t^-i ^ 

It follows from r — / £ N that for each 77 = 0, . . . , [ZJ we have r — Z + 77 G N and 
Q<r — l + 'q<r — I, which implies that 

(15) \u\r-i+n..p.n<\Mr,p,n- 

Combining the preceding bounds we obtain for all d < 0^ and = 0, . . . , [ZJ , Z that 



(16) 



n 9^ 

iai<M 



If g = cx) then the result follows immediately. If 1 < g < 00 it follows from JT 
that 

\W\\i,,,n < C-K-i[l]+2f'^[d^-'-/P+-/'>\\u\l^^. 

_j_^"/7+A'-' 



n ^""1 
iai<M 

The result then follows by incorporating the constant +2 
dependence on [l\ substituted with n,r,p and q. 

Only the case that p = q = m: ^ 2 will be used in Section [S] 



into C, with the 
□ 



Corollary 2.10. Let r be a real number, ^ be a nonnegative integer such that 
r — fi > n/2. Then, there exist three positive constants 0^ (dependent on 9,p,n 
and r), C (dependent on ^l,n and r), and K (explicitly dependent on \r\ — I, 
cf. pU]) with q = 2) satisfying the following property: for any set A C ^l, such that 
d — S [A, fl) <dr, u e W^'"^ (ri) and real number I £ [0, r — ^] such that r — / e N, 



(17) 



jr — l 



n 9^ 

\a\<fi 



3. Application: Unsymmetric Meshless Methods for Operator 

Equations 

In this section, only the dependence of constants on the discretization parameters 
r and s is explicitly stated since we regard the spaces and mappings involved as 
fixed. 

We now apply the sampling inequality stated in Corollary 12.101 to Schaback's 
framework for unsymmetric meshless methods for operator equations |12| . Due to 
unaddressed issues contained in its original formulation, we use a modified version 
stated in this section. On a technical level, it differs substantially, e.g., certain 
spaces have been eliminated, the inequalities apply over possibly different spaces, 
and the proof of the error bound has been slightly modified, but the underlying 
ideas are the same and entirely due to Schaback [12] . 
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The framework provides an error bound for meshless methods which approxi- 
mately solve a linear operator equation in the following setting. The first require- 
ment is a continuous and bijective linear operator L : U F mapping from the 
solution space to the data space. The spaces U and F are assumed to be complete 
in order to ensure the boundedness of L^^ : F ^ U. It is also assumed that the 
exact solution u* d U where U C U is called the regularity subspace. We denote 
F := LU. The framework requires a scale of finite-dimensional trial subspaces 
Ur a U equipped with a projector Ilr : U Ur- The framework requires a lin- 
ear, continuous, and bijective test mapping A : — > T, where the test space T is 
assumed to be complete in order to ensure the boundedness of A~^. We denote 
T := AF. Test data from T is discretized into finite-dimensional test subspaces Ts 
with a test discretization mapping 



(18) 



IT, : T ^ 



the operator norm of which must be bounded by a constant, which is independent 
of s. It follows that the operator norm of 



7r,AL -.U^n 



is bounded similarly since 

ll7r,AL 



< IItt, 



l|AL||, 



In order to apply the error bound of Schaback's framework a number of inequal- 
ities must be supplied. The first of these is the trial space approximation property 



(19) 



\u — Ilru\\^ < e (r) \\u\\jj for all u e U. 



The second inequality is the test discretization's stability condition 



(20) 



1" 1 1 T. 



for all Ur G Ur 



If the stability factor /? (s) grows as the test discretization is refined, i.e., as s 
decreases towards zero, then the order of convergence in the final error bound (j23p 
will be less than that provided by the trial space approximation property (|19p . 
When the stability factor does not grow, the test discretization is called uniformly 
stable. The final inequality required by Schaback's framework involves a numerical 
method capable of providing an approximate solution u* ^ G Ur which satisfies the 
numerical method approximation property 



(21) 



TTcA (Lul 



f)\\^ <C||7r,Ai||^_ 



In particular, if the numerical method computes u* ^ 
hand side of ((2T|) then the constant is at most one, since 



e(0 Ik* lie/ • 

G Ur which minimizes the left 



|7r,A(Lu;, -/)| 



(22) 



< 
< 
< 



llvr.Ai {nru*-u*)\\, 



\nsAL\ 
k.Ail 



u- 



.T,<^ir) II"* lb- 



Theorem 3.1. [12' Theorem 1] Given the setting stated above, if the inequali- 
ties (|19|) . ()20p . and (|2ip are satisfied then the following error bound holds: 



r,s\\ij 



< 



1+Pis) (AL) 



k.^Ai 



il + C))eir)\\u*\ 
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Proof. We have that 



UrU* — U* ,, 

< eWI|w*|ly+||n,M*-<J^ 



|n,.M* — u*.s||^ < 



(AL) ^ \\AL {nru* ~ u;,)\ 



TT.AL (11,. 



T— (7 



r.s y 1 1 2^ 



TTgAL (u* — u. 



r.s } 1 1 7^ 



Hs) (AL)-i ^^Jn,AL\\^^^^e{r)\\u*\\^{l + C) 



< ^{r)\\u*\\u 
+/3(s) (AL) 



l7r,AL||^_^j,^6(r)|i^.*|ij^(l + C), 



□ 



The stabihty condition pop can be estabUshed using an inverse estimate 

(23) \\ur\\tj < 1 {r) \\ur\\jj for all Ur £ C/^, 
a sampling inequality 

(24) II/IIj, < C (a (s) ||/||^ + /3 {s) |k./||^J for all / G f , 
and ensuring that a /ine enough test discretization is chosen such that 

(25) " - ■ ' ' 



Ca{s)j{r)\\AL\\^^f {ALy' 



T^u - 2' 



where C is the constant appearing in (|24p . Typically, 7 (r) ^ oo as r ^ 0, while 
a (s) ^ as s ^ 0. 



Proposition 3.2. [12, Theorem 2] // and ^ hold then so does (|20l). 

Proof. We have that 

||ALu,.||y < C (a{s)\\ALur\\f+ P{s)\\TrsALur\\rpJ 

< C{a (s) ||AL||^_^ \\ur\\o + P (s) ||7r,Aiu,|| ^.J 

< C{a (s) \\AL\\^_f 7 (r) \\ur\\u + (3 {s) ||7r,AL^.,||^J 

< c(a{s)\\AL\\^^^^{r) (Ai)"^ \\ALu4^ 

+P{s) Wir^ALurW^J 

< ^\\ALUr\\T + Cf3{s)\\TTsALUr\\rr^ 

and the result follows by incorporating the constant 2C in /? (s). □ 
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3.1. Convergence Results for the Poisson Problem. We consider the example 
from [12|, Section 4.1], a Poisson problem with mixed, inhomogeneous boundary 
data: let f2 be a bounded domain in with a Lipschitz-contimious boundary. We 
denote fli := fl, fl2 ■= C dil, and fl^ = C dil so that the dimension of 
each domain is given by ni = n, and n2,n3 — n — 1. Let m, rh be nonnegative real 
numbers such that to — m e N, and 



(toi,TO2,TO3) 

U 
F 



= (to, TO + 3/2, m+ 1/2) 



(26) Lu := ( —Au.ulrD , -^\yn 

\ on 

with analogous definitions made for (toi,TO2,TO3), U, and F. With the space F 
equipped with the norm ||-||^ :— \\-\\pi + \\-\\p2 + \\-\\p3, it follows that the hnear 
operator L, as defined above, is continuously invertible either as L : U — > F or 
L:U ^F. 

We assume that the solution comes from U and that the trial space Ur is chosen 
such that the trial space approximation property holds with e{r) = O (r™"™) , 
a property satisfied by kernel-based meshless trial spaces, c.f. Narcowich et al. [3 [5], 
and finite-element trial spaces [3, Theorem 4.5.11]. We also assume that the inverse 
estimate (|23l) holds with 7 (r) — O (r™~™), as is the case for finite-element trial 
spaces [31 Theorem 4.4.20]. Obtaining an inverse estimate with the expected factor 
7 (r) =0 (^|.™-™^ appears to be an open problem for kernel-based meshless trial 
spaces. Narcowich et al. provide an inverse estimate with the expected factor 
for the case of Sobolev spaces over R". Both Schaback and Wendland ITT], and 
Duan [4 provide inverse estimates for Sobolev spaces over a domain. Unfortunately, 
the factor involved in these inverse estimates are worse than the finite-element 
case. Further progress on this problem is expected to be reported in the thesis of 
Rieger 0. 

We consider the case of strong testing here, which means that the test mapping 
A : F ^ T is just the identity mapping and that each test space T'^ coincides with 
the corresponding data space F'^. Weak testing is also possible, in which case the 
test functionals integrate functions in F'^ against test functions, resulting in the 
test data in each T*^ acquiring additional smoothness. This is discussed in detail 
by Schaback [121 [14]. Each domain is discretized onto finite subsets C ftk, with 
the same fill distance s = S (Y^, 57^). Furthermore, they are assumed to satisfy the 
property that is bounded by s^"* up to a constant, as is the case for domain 
discretization with a uniformly bounded mesh ratio [T^]. We define discrete test 
spaces 

{Q:|a|<Mfc}-#>;'' 



equipped with a norm 
(27) 

and a test discretization tt^ -.T^ ^ 

(28) ^Ifk := n d'^hW^ for aU A e T 



k 
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where /ife is an integer such that mk — fJ.k — nk/2 > 0, and furthermore this difference 
is independent of k. The discrete test space Ts :— x x is defined and 
equipped with a norm, analogously to F and T. The test space T is then equipped 
with a test discretization 1:^ : T —>■ Tg defined by 

(29) ttJ (^1/1,^,^2, TT^/s) for aU / = (A, A, /a) e T, 

Proposition 3.3. If for each k, rrik—fJ-k—nk/'^ > then tTs : T Tg is well-defined 
and the operator norm \\t^s\\t^t bounded independently of s. 

Proof. Suppose / = (/i, 727/3) £ T. Since ruk — fik > nk/2 we have from the 
Sobolev embedding theorem that T*^ ^ C^* (Ofc) and therefore the test dis- 
cretization is both well-defined and there exists some constant independent of 
/ = (/i, /2, /a) such that for each fk, 

ll/fc|lcffc(?2^) - C'll/fcllT" • 

Since is bounded by s~"'= up to some constant which is independent of s, it 

follows that 

k=l 

= E^"'' E E 

k=l x£Y^ |Q|<^fc 
3 

\ ^ ^rii. II ^112 



< E ll/llc^^ (nl) # {« : l«l < l^k} ■ 



< c-^ 11/11?,. = C 11/11?,, 
fc=i 



□ 



Proposition 3.4. There exists a constant sq such that for all s < Sq a sampling 
inequality (|24|) /lo/rfs iwit/i a constant C for the test space T and test discretization 
TTs-.T with a (s) s""™, and [i (s) := s/^i-™i = s;'2-i/2-m2 = gM3-i/2-m3 _ 

Proof. From Corollarv l2 . 1 01 and (^7)) . it follows that for each k there exist constants 
Ck and Sfe such that for s < sq min (1, si, S2, S3), 

IIMIt^ < Ck[ais)\\fk\\f.+s^--"^-\\n':fk\\^,) 
< Ck(ais)\\fk\\f,+s^--"^^\\n':f4^,) 

since fii — mi < M2 ~ 1^2 = Ms ~ '^3- The result then follows with a constant C by 
combining the preceding inequalities. □ 

We assume that the s is sufficiently small to satisfy the requirements of Propo- 
sition and (pS)) . Even in the fractional case, the sampling inequality introduced 
here provides a (s) = g™""* which shrinks as rapidly as the expected inverse esti- 
mate factor 7 (r) = r™~™ grows and thus s and r can be kept proportional. This 
is in contrast to previous sampling inequalities which necessarily introduce a fac- 
tor a (s) = when bounding fractional order Sobolev norms, requiring the 
test discretization to be refined more rapidly than the trial discretization and thus 
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Ml =0 


None 


m — m — 2 


m — m — 3 


rn — m — 4 


Ml = 1 


None 


None 


rh — m — 2 


rh — m — 3 


Ml -2 


None 


None 


None 


TO — TO — 2 


Ml = 3 


None 


None 


None 


None 



Table 1. Order of convergence in various Sobolev norms estab- 
lished by a modified formulation of Schaback's framework, using 
trial spaces with optimal properties and strong testing with var- 
ious order test discretizations to solve two- or three-dimensional 
Poisson problems. 



diminishing the order of convergence by [to] — to. If the function u* ^ G Ur which 
minimizes the left hand side of pT|) has been computed, then Schaback's frame- 
work provides the error bound with constant C — 1. The order of convergence 
established by this error bound, in terms of both the trial and test discretization, 
is then given by (3 (h) e (h) and satisfies 

Table [T] states particular convergence results using various Sobolev norms and 
test discretizations in two- and three-dimensions. These results show that, for 
convergence in higher order norms, the highest order of convergence is obtained 
using a higher order test discretization introduced here. 

We note that to obtain a uniformly stable test discretization with /? (s) = 1 
would require choosing the order Mfc of each test discretization to be equal to that 
of the order mk of the test space. Unfortunately, this does not seem to be possible: 
in order for the test discretizations operator norm to be bounded independently of 
s, the order of each test space is required to be greater than that of the test 
discretization Mfe by at least nj,/2. It follows that the order of convergence, in terms 
of both the trial and test discretization, provided by this modified formulation of 
Schaback's framework is always less than that of the trial space approximation 
property. Another consequence is that the order of U must be at least 2 -I- n/2, 
and therefore convergence in the norm can only be concluded suboptimally 
from convergence results in higher order Sobolev norms, using strong testing in 
this modified formulation of Schaback's framework. 

4. Conclusions 

We have further generalized the sampling inequalities of Arcangeli et al. [T], 
Madych [6], and Wendland and Rieger T5], to optimally bound fractional order 
Sobolev semi-norms, and to incorporate higher order data into the bound. When 
used in a modified formulation of Schaback's framework to prove convergence rates 
for unsymmetric meshless methods this new sampling inequality has two benefits: 

(1) It results in more optimal estimates for problems involving fractional order 
Sobolev spaces, particularly by providing a more optimal constant a (s). 

(2) For convergence in higher order Sobolev norms, higher order results are 
obtained using a higher order test discretization in comparison to the zero 
order test discretization. 
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The zero order test discretization has been widely employed in practice, and corre- 
sponds to what is usually called Kansa's method or unsymmetric collocation. On 
the other hand higher order testing has not, and its value in practical applications 
requiring convergence in stronger norms is an open question worthy of further study. 
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